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Resumo 



Nesta tese sao estudados varios proccssos que envolvem buracos negros em quatro e mais dimensoes. 
Primeiro, usando tecnicas de teoria de perturbagoes, compara-se a radiagao escalar sem massa e com 
massa emitida por uma particula altamente energetica caindo radialmente dentro de urn buraco negro 
de Schwarzchild. Mostra-se que em tempos muito grandes, o sinal de perturbagoes escalares massivas 
e dominado por uma cauda oscilatoria universal, que aparece devido a efeitos de curvatura. Mostra-se 
tambem que o espetro de energia esta de acordo com o calculo do ZFL depois de considerar o teorema 
do nao-cabelo. Na scgunda parte, estuda-se o fcnomeno de superradiancia em dimensoes genericas e 
conjectura-se que a energia maxima extrai'da de um buraco negro em rotagao pode ser compreendida 
em termos do volume proprio da ergoregiao. Finalmente, estudam-se algumas consequencias do efeito de 
superradiancia na dinamica de luas que orbitcm em torno de buracos negros com rotagao em ccnarios com 
mais de quatro dimensoes. Em quatro dimensoes, luas em torno de buracos negros geram mares de baixa 
amplitude, e a energia extrai'da da rotagao do buraco negro e sempre menor do que a radiagao gravitational 
perdida para o infinito. Mostra-se que, em dimensoes maiores do que cinco, a energia extraida do buraco 
negro atraves da superradiancia e maior do que a energia emitida para o infinito. Estes resultados dao 
um forte apoio a conjectura de que cstc cfcito e a rcgra, e nao a excepgao, em dimensoes superiores. A 
superradiancia domina e luas espiralam para fora; para uma dcterminada frcquencia orbital, a energia 
extrai'da no horizonte e igual a energia emitida para o infinito e " orbitas flutuantes" ocorrem. Da-se uma 
intcrpretagao desse fenomeno em termos do paradigma da membrana e da aceleragao das mares devido 
a dissipagao de energia no horizonte. 

Parte dos resultados obtidos durante esta tese figuram na Ref. Q]. 



Pa.la.VraS-Cha.Ve; Buracos negros; Campos escalares; Superradiancia; Efeitos de Mare; Di- 
mensoes extra. 
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Abstract 



In this thesis we study several dynamical processes involving black holes in four and higher dimensions. 
First, using pcrturbative techniques, we compare the massless and massive scalar radiation emitted by a 
particle radially infalling into a Schwarzchild black hole. We show that the late-time waveform of massive 
scalar perturbations is dominated by a universal oscillatory decaying tail, which appears due to curvature 
effects. We also show that the energy spectrum is in perfect agreement with a ZFL calculation once no- 
hair properties of black holes are taken into account. In the second part, we study the phenomenon of 
supcrradiance in higher dimensions and conjecture that the maximum energy extracted from a rotating 
black hole can be understood in terms of the ergoregion proper volume. We then study some consequences 
of supcrradiance in the dynamics of moons orbiting around higher-dimensional rotating black holes. In 
four-dimensional spacetime, moons around black holes generate low-amplitude tides, and the energy 
extracted from the hole's rotation is always smaller than the gravitational radiation lost to infinity. We 
show that in dimensions larger than five the energy extracted from the black hole through supcrradiance 
is larger than the energy carried out to infinity. Our results lend strong support to the conjecture that 
tidal acceleration is the rule, rather than the exception, in higher dimensions. Superradiance dominates 
the energy budget and moons "outspiral" ; for some particular orbital frequency, the energy extracted 
at the horizon equals the energy emitted to infinity and "floating orbits" gencrically occur. We give an 
interpretation of this phenomenon in terms of the membrane paradigm and of tidal acceleration due to 
energy dissipation across the horizon. 

Part of the results obtained during this thesis appeared in Ref. [lj]. 



Keywords; Black holes; Scalar fields; Superradiance; Tidal effects; Extra dimensions. 
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Chapter 1 



Introduction 



This thesis is devoted to the study of several processes involving black holes in four and higher dimensions. 
In particular, we shall discuss processes that involve the emission of scalar radiation by black holes in 
four and higher dimensions within a perturbative approach, i.e. by solving the linearized field equations. 
The first part of the thesis will be devoted to the signatures of high-energy processes around non-rotating 
black holes, when ordinary matter is coupled to a massive scalar field. In the second part, we shall discuss 
strong tidal effects which are related to the phenomenon of superradiance in higher dimensional black 
hole spacetimes. 

Black holes are the simplest macroscopic objects in the Universe. They are easily understood using 
only the concept of space and time given by General Relativity. Furthermore, they can be described 
by a couple of parameters, namely, their mass, angular momentum, and charge. This is the famous 
no-hair theorem, which states that, in its final state, a black hole is uniquely described by these three 
parameters 4]. In the last fifty years there have been considerable progresses in the understanding of 
black holes, and the mathematical tools to describe them have been highly improved [3]. Nowadays, it 
is universally recognized that black holes are not only of academic interest, but they are also of central 
importance in astrophysical processes and fundamental physics. For example, it is believed that most 
galaxies contain supermassive black holes at their centre [J] and that the formation of astrophysical black 
holes is probably related to extreme phenomena. In fundamental physics they play a key role since 
Hawking's semi-classical prediction of black hole evaporation by emission of a thermal radiation due to 
quantum effects Hawking's radiation was the first phenomenon to be predicted considering both 
gravitational and quantum effects. Black holes may therefore play a major role in the attempt to find a 
consistent quantum theory of gravity. 

More recently, the discovery of the gauge/gravity duality has given a whole new interest to the 
study of general spacetimes, and more particularly of black hole spacetimes. This duality maps the 
dynamics of some strongly coupled quantum field theories (QFTs) in D-dimcnsional spacetimes to the 
dynamics of semiclassical gravity in D + 1-dimensional spacetimes. It is actually the only tool available 
to study the dynamics in strongly coupled QFT The most widely studied gauge/gravity duality 
is the AdS /CFT correspondence, which maps the dynamics of non-Abelian conformal field theories 
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(CFTs) in D-dimensional spacetimes onto semiclassical gravity in asymptotically Anti-de Sitter (AdS) 
spacetimes, i.e. , spacetimes with a negative cosmological constant, in D + 1 dimensions [3, H|- The field 
theories described by this correspondence are very different from the most important non-Abelian theory 
describing the physical world, Quantum Chromodynamics (QCD), the theory that describes the strong 
nuclear force. However, there is great hope that there could be a similar correspondence for some features 
of QCD (e.g. [9-12|)_ and some recent results obtained at the Rclativistic Heavy Ion Collider (RHIC) 



reinforced this hope 



13j . This duality could provide tools to understand the behavior of condensed matter, 



quark-gluon plasmas, and other strongly coupled systems, impossible to describe using perturbativc 
methods, in terms of AdS black holes that interact semiclassically with fields. 

In addition to the AdS /CFT duality, in the last decades, there has been a growing interest in physical 
phenomena in higher dimensional spacetimes, mainly motivated by higher dimensional solutions which 
naturally arise in the context of string theories and supergravities. In some scenarios, the extra dimensions 
arise naturally as an attempt to solve the hierarchy problem. Put in simple words, this problem seeks to 
understand why there is such a big difference between the electroweak scale tuew ~ 300 GeV and the 
gravity or Planck scale M^pi ~ 10 19 GeV. In some extra-dimensional models, if the extra dimensions are 
highly warped and correspond to a very large volume, the Planck scale can be as low as Mdpi ~ 0(1) 
TeV, eliminating the large difference between this scale and the electroweak scale. These are the so-called 



TeV-gravity scenarios 



14 



17j . The idea behind that comes from string theory and the concept of p- 



brancs. A p-brane is a hypersurfacc with p-spatial dimensions. Gravity would live in a D-dimensional 
spacetimc, and the gauge interactions, i.e. , the electromagnetic and nuclear forces, would be constrained 
to a 3 + 1-dimensional brane. The gravitational attracti ve p ower would then be "diluted" in the extra 
dimensions, appearing to be much smaller in the brane fol. 1 1 sjj - Presently the LHC has begun to search for 
evidences of Tev-gravity models, looking for trans-Planckian signatures, i.e. , processes where the energy 
involved exceed the Planck energy. At the Planck scale, non-linear quantum gravity effects dominate and 
we would need a complete quantum gravity theory to describe the physics near this domain. However, in 
some conditions gravity becomes highly non-linear, and black hole formation is expected 19M22I]. This 



happens when the impact parameter b = J/y/s, with J being the angular momentum and y/s the center- 
of-mass energy, is of the order of the Schwarzschild-Tangherlini radius R s . which is a generalization of 
the 4-dimcnsional Schwarzschild radius R s = 2M in higher-dimensions. The formation of these higher- 
dimensional black holes would carry a clear signature through the decay by Hawking radiation, emission of 
gravitational waves, and in some theories, other types of radiation, such as scalar radiation. Although at 
the time of writing, none event of black hole formation has been recorded at the LHC 23|,|24| , further study 



is necessary to fully understand the trans-Planckian regime and the signatures that higher-dimensional 
black hole formation could carry. 

Today it is widely accepted that there are at least three complementary approaches to solve the Ein- 
stein's equations for a given problem: (i) analytical computations, (ii) semianalytical methods and soft 
numerics, (iii) fully numerical simulations. Purely analytical computations are normally only possible 
at a linear level and, in most cases, it is impossible to obtain a full analytical solution. In fact, solving 
the Einstein's equations exactly, i.e. , with no approximations, is a formidable task which requires super- 
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computers and sophisticated numerical methods. Complementary to these, semianalytical methods are 
instrumental for a better understanding of dynamical processes in General Relativity and sometimes for 
interpreting the results of the simulations. There is a need for the development of new approximation 
schemes to accompany large scale simulations. Fully numerical simulations give only numerical answers 
to problems. Often these answers do not provide physical understanding, or even what principles are 
behind the process. Furthermore, in some cases, simulation results can be simply incorrect or misleading. 
By closely coupling various perturbation schemes it is possible to interpret and confirm simulation results. 
On the other hand, even using a perturbative approach usually requires numerical solutions. The devel- 
opment of these complementary approaches cannot be done independently. This is the main motivation 
of this thesis, to study and understand processes that could be seen in future numerical simulations and 
could also be of great importance, both in the astrophysical context and in more theoretical applications, 
like the gauge/gravity duality and TeV-gravity scenarios. 



1.1 Outline of the thesis 

This thesis is organized as follows. In chapter [5] we compute the fundamental equations of scalar radia- 
tion in terms of black hole perturbations sourced by a test-particle in geodesic motion around a spinning 
black hole. We derive the Teukolsky equations and solve the wave equation using the Green's func- 
tions approach. In chapter [3] we specialize the problem to the case of a particle falling radially into a 
Schwarzschild black hole. We compute the massive scalar radiation waveforms emited by the particle and 
the correspondent energy spectra. In chapter [4] we study the superradiant scattering of a massless scalar 
field with a singly spinning black hole in D = 4 + n dimensions. In chapter [SJ extending the discussion of 
Ref. [25)], we study the effect of an electrically charged particle orbiting a neutral central object in 4 + n 
dimensions. By applying the membrane paradigm |26[, we derive a simple formula for the ratio between 
the energy flux at infinity and at the horizon in 4 + n dimensions. Furthermore, we confirm the results 
obtained at the Newtonian level by solving the wave equations using black hole perturbation theory both 
analytically and numerically. We conclude in chapter [5] 

Throughout the thesis we use G = c = h = 1 units, except in the first section of chapter [5] where, for 
clarity, we show G and c explicitly. 
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Chapter 2 



Scalar perturbations of singly 
spinning Myers— Perry black holes 



The way a black hole reacts to external perturbations provides us with a deep understanding of the 
space-time around it. In this thesis we will be interested, more particularly, in scalar perturbations of 
singly-spinning Myers-Perry black holes due to the presence of a test particle coupled to the scalar field. 

Scalar fields are predi cted b y several theories aiming to unify Einstein's theory with the Standard 
Model of particle physics 



27 



- 30| ■ Their existence is being widely looked, and therefore, the study of such 
scalar fields is very relevant and timely. 

We will not go into details about the mathematical formalism of the problem, since it as been deeply 
studied over the years and it is not the purpose of this thesis (for further detail about black hole pertur- 
bation theory see Q). The starting point of a rclativistic perturbation theory is to consider two different 
spacetimes, the physical one, which carries the information about the actual physical system, and the 
background one, which corresponds to a simpler idealized spacetime, solution of Einstein's equations. The 
background and physical spacetimes can be related using different maps, each of them corresponding to 
different a gauge choice. The transformation between maps is then called gauge transformation. For a 
particular map, we can write the spacetime metric as 

g ab (xn=9^(xn + h ab (x»), (2.1) 

where (a; M ) is the background metric, and h a b{x^) is a small perturbation that can be obtained solving 
Einstein's equations order by order in a small curvature and small velocity expansion. To avoid possible 
physical interpretation problems due to the choice of the gauge, one usually looks for gauge-invariant 
quantities, i.e. , quantities that have the same value independently of the chosen gauge. 

First-order perturbations of black holes began with the pioneering work of Regge and Wheeler 3l| , 



who studied the stability of Schwarzschild black holes, followed by the work of Zerilli, who first gave 
a fully rclativistic treatment of the gravitational radiation emitted by a particle falling radially into 



a Schwarzschild black hole 



Finding a couple of decoupled master equations describing generic 
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perturbations in the case of a rotating black hole was much more difficult to solve, but scalar perturbations 
were shown to be easily treated even in this case jsjj, [3] . The separability of the scalar wave equation 
in the Kerr family of geometries, as demonstrated by Carter in 1968 [35 1, and explicitly exposed by Brill 
et al in 1971 j^ . opened a new door for the study of rotating black hole perturbations. Among many 
others, one of the major achievements of this era, the so-called "golden age" of General Relativity, was the 
discovery, by Tcukolsky in 1973 [36| . of fundamental equations governing, not only scalar perturbations, 
but also, gravitational, electromagnetic, and massless fermionic perturbations. 

In this chapter we shall generalize some of these results to higher-dimensional backgrounds, considering 
a test-particle coupled to a massive scalar field around the spinning black hole. The results here derived 
will serve as a framework for the subsequent chapters. 



2.1 The background metric 



In four dimensions, there is only one possible angular momentum parameter for an axisymmetric space- 
time, and rotating black hole solutions are uniquely described by the Kerr family. In higher dimensions 
there are several choices of rotation axis, which correspond to a multitude of angular momentum param- 
eters 37J . Here we shall focus on the simplest case, where there is only a single axis of rotation. In the 



following we shall adopt the notation used in Rcfs. 



4C|, to which we refer for details. 



The metric of a (4 + n)-dimcnsional Kerr-Myers-Perry black hole with only one nonzero angular 



momentum parameter is given in Boycr-Lindquist coordinates by 
ds 2 = 



37] 



A - a 2 sin 2 1? , , 2a(r 2 + a 2 - A) sin 2 , , J (r 2 + a 2 ) 2 - Aa 2 sin 2 . , n , 
dt 2 dtdep + - '—- sm 2 -ddcj) 2 



A 



dr 2 + Sdtf 2 + r 2 cos 2 tidSt, 



(2.2) 



where 



r + a cos ■& , 



2Mr 



l-n 



(2.3) 



and dfl 2 t denotes the standard line element of the unit n-sphere. This metric describes a rotating black hole 
in asymptotically flat, vacuum spacetime, whose physical mass M. and angular momentum J (transverse 
to the rtfi plane), respectively read 



M= {n + 2)An+2 M. 



8- 



J = 



-Ma . 



(2.4) 



where A n+2 = 27r" +3 / 2 /r[(n + 3)/2]. 



The event horizon is located at r = rn , defined as the largest real root of A. In four dimensions, an 
event horizon exists only for a < M. In five dimensions, an event horizon exists only for a < s/2M, and 
the black hole area shrinks to zero in the extremal limit a — » \2M. On the other hand, when D > 5, 
there is no upper bound on the black hole spin and a horizon exists for any a. 
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2.2 The wave equation 

We consider a small object in a geodesical curve around a spinning black hole and a scalar field of mass 
rn s = fih coupled to matter (from now on we set h= 1. In these units p, has the dimensions of 1 /length). 
At first order in perturbation theory, the scalar field equation in the background (|2.2[) reads 



(2.5) 



where a is some coupling constant. For simplicity we focus on source terms of the form 

(It 



T=- 



-9 



5 (i+n \x - X{t)) , 



(2.6) 



which corresponds to the trace of the stress-energy tensor of a point particle with scalar charge q p . 

Because of the coupling to matter, the object emits scalar radiation, which is governed by Eq. (|2.5 
To separate Eq. (|2.5[) . we consider the ansatz 



<p(t,r,4,<f>) = J due* m +- iut R(r)S lmj (#)Yj, 

l,m,j 



(2.7) 



where Yj are hyperspherical harmonics 
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41| on the n-sphere with eigenvalues given by + n — 1) 



and j being a non- negative integer. The radial and angular equations read 



-n d ( ,i a dR 
r n — r n A — 
dr V dr 



[uj(r 2 + a 2 ) - ma] 2 j(j + n— l)a 2 



A - [i 2 r 2 \ R = T lmj 



and 



1 



sin ■& cos™ ■§ dfl 



d I . „ a dSi m j 

sin v cos v - 

dv 



{uj 2 - [i 2 )a 2 cos 2 z? 



i(j + n-l) 



sin 2 d cos 2 t? 



-4, 



Imj 



Slmj — . 



(2.9) 



where A = Ai m j — 2mu>a + u> a , A\ m j are the eigenvalues of the angular equation, and we have defined 



l,m,j 



Radial equation 

Defining a new radial function Xi m j (r) 



(2-10) 



X lmj = r n / 2 (r 2 + a 2 )^ 2 R, 



(2.11) 



we get the nonhomogencous equation for the scalar field 



dr 2 



V 



{r 2 + a 2 ) 3 / 2 



(2.12) 



where drjdr* = A/(r 2 + a 2 ) defines the standard tortoise coordinates and the effective potential V reads 



7 



V =uj 2 



3r 2 A 2 A[3r 2 +a 2 -2Mr 1 -"(2-n)] 



1 



r 2 + a 2 (r 2 + a 2 ) 4 



(r 2 + a 2 ) 3 



(r 2 + a 2 )' 



{a 2 m 2 - AM 



amuj . / 9 9 , \ 



n(2-n)A 2n(l-ra)M j(j + ?i-l)a 2 

77- ~t 



4r 2 



2r™ +1 



(2.13) 



Angular equation 

In the low-frequency limit the angular Eq. (|2.9[) can be solved exactly. For the massless case /j, = and 



at first order in auj, the eigenvalues can be computed analytically 



41] 



A kjm = (2fc + j + |m|)(2fe + j + \m\ + n + 1) + 0(au) 



(2.14) 



By setting 2k = I — (j + \m\), the eigenvalues above take the form Aij m = 1(1 + n + 1) and I is such 
that I > (J + \m\), which generalizes the four-dimensional case. An important difference from the four- 
dimensional case is that regularity of the angular eigenfunctions requires k to be a non-negative integer; 
i.e. for given j and m only specific values of I are admissible. In fact, it is convenient to label the 
eigenfunctions and the eigenvalues with the "quantum numbers" (k,j,m) rather than with (l,j,m) as in 
the four-dimensional case. Note that when aio > 1, these eigenvalues might not be accurate. Therefore, 
in some cases, we will use the exact numerical eigenvalues. These can be obtained assuming an expansion 



for the eigenfunctions of the form 



4l|, 



S kjm = sin(tf)l m l cos(t?) J ' ^ 5 p (cos 2 . 

p=0 



(2.15) 



This series (if convergent) automatically satisfies the regularity conditions at 8 = 0, tt/2, it. Upon substi- 
tution of (|2.15p into the angular equation (|2.9|) . we obtain a three-term recursion relation 



a di + f3 a = , 

a p a p+ i + f3 p a p + >a p _i = , (p = 1, 2, . . .) 



(2.16) 



where 



a p = -2(p + l)(2j + n + 2p+l), 

h = U + H + 2 P)U + n + H + 2p + 1) - A kjm , 

7 P = -(auj) 2 . 



(2.17) 



Then, given a value for aw, the eigenvalues Akj m can be obtained solving numerically the continued 



fraction equation 



42] 



/3o 



"071 



a 172 



= 0. 



(2.18) 



The (non-normalized) zeroth-order eigenfunctions are given in terms of hypergeometric functions 
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|4l| 



S kjm oc sin{^ ml x J F -k,k + j + \m\ + + ; x 2 



(2.19) 



where x = cos($). We adopt the following normalization condition: 

f.?r/2 

did sin $ cos $Skj 7n S , 



kjm J kjm 



1 



(2.20) 



where the integration domain has been chosen in order to have a nonvanishing measure also in the case 
of odd dimensions. Note that this normalization differs from that adopted in Ref. 41 [. 



2.3 Green function approach 

To solve the wave equation, let us choose two independent solutions X r k " m and X£° m of the homogeneous 
equation, which satisfy the following boundary conditions: 



yoo 
kjm 

V r H 

^kjm 



r — > oo 



out 6 



> r >+A h 



(2.21) 



Y"oo 
kjm 

X™ ~ e 



outt 

— iknT^ 



(2.22) 



Here fc# = w — rniln, fcoo = \/u 2 — /i 2 , and = — lim,._>. rff gt4>/g<p<t, = a/{i" 2 H + a 2 ) is the angular 
velocity at the horizon of locally nonrotating observers. The Wronskian of the two linearly independent 
solutions reads 



W = X 



r H aA kjm 



-XI 



dX kjm 



(2.23) 



Lfe j' m dr, " kjm dr, 

and it is constant by virtue of the field equations. 

Imposing the usual boundary conditions (see Appendix |A"]) . i.e. , only ingoing waves at the horizon 
and outgoing waves at infinity, Eq. (|2.12[) can be solved in terms of the Green function 43 1 



X kjm (r,) 
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kjm 
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A r 'n/2 

J 1 ,. ( r ') — X th dr' 

±kjm,V ) , 2 2)3/2 ^kjm ar * 



^kjm 

w 



Ar /n/2 



r . Tk]m{r '\ r ' 2 + a 2)3/2 X S 



X°° dr' 



(2.24) 



In the next chapters we shall use these results in two different scenarios: a particle radially infalling 
into a Schwarzschild black hole in four dimensions, where the particle is coupled to a massive scalar field; 
and a particle coupled to a massless scalar field in an equatorial circular geodesic around a singly spinning 
Myers-Perry black hole. 
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Chapter 3 



Scalar radiation from an infall of a 
particle into a Schwarzschild black 
hole 



In this chapter we will use the results presented in the last chapter in a particular scenario, a particle 
radially infalling into a Schwarzschild black hole coupled to a massive scalar held. All the formulae 
derived in the last chapter apply, as long as we use a = and n = 0, i.e. , a non-rotating black hole in 
the usual four dimensions. 

In the context of TeV-scale gravity and gauge/gravity duality scenarios, the signature of these kind 
of collisions could be extremely relevant. Although we restrict to the non-rotating case and to four 
dimensions, it would be very interesting to generalize the results derived below to rotating black holes in 
higher dimensions, and even to non-asymptotically flat geometries. 

The use of perturbative techniques to study the energy radiated by an infalling particle into a 
Schwarzschild black hole goes back to 1970 with the work of Zcrilli [32 1 and Davis et al 



44 



who first 



computed the gravitational energy radiated away by a small test particle of mass m p , fal ling radially 



45| has con- 



from rest at infinity, into a Schwarzschild black hole of mass M. A recent work by Mitsou 
firmed and improved significantly the numerical accuracy of the results obtained by Davis et al. Later, 
Ruffini generalized these results, allowing the particle to fall with an initial velocity at infinity 46|. The 



importance of the infalling particle model lies in the fact that it sometimes appears as a limit case of 
more general scenarios such as the coalescence of black hole binaries in the extreme mass-ratio limit 47 1 . 



Furthermore, the limit m p — > M describing the collision of two black holes, do predict reasonable results 
still with in p erturbation theory, making perturbation theory a fundamental tool to study important phc- 



47 



48| . Recently, all these results where extended to the case of large boost factors by Cardoso 
and Lemos (4J, where they considered a massless particle falling radially from infinity. 

One interesting common behavior in this kind of processes, in the case of massless fields, is that the 
black hole fundamental quasinormal frequency acts as a cutoff in the energy spectra. Furthermore, the 
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signal is always dominated by the quasinormal ringing at intermediate times. Quasinormal modes are the 
characteristic oscillations of a black hole. They are completely independent of the initial configuration 
that caused such vibrations. Therefore, they are a characteristic of black holes. The name derive from 
the similarity between these excitations and normal mode systems, such as the normal modes of a guitar 
string. However, they are called quasinormal for two reasons: first, they are not stationary modes, since 
they are exponentially damped; secondly, unlike a normal mode system, quasinormal modes seem to 
appear only over a limited time interval, at very late times the quasinormal ringing gives way to a power- 
law falloff. Mathematically, this is related to their incompleteness. For full reviews about quasinormal 
modes in different black hole spacetimes see for example Refs. 50l452|. 



Here we will be interested in the energy radiated due to the presence of a massive scalar field. Massive 
scalar fields present a number of interesting behavior in the presence of black holes. For instance, it is 
known that the quasinormal decay of massive scalar fields is slower than massless scalar fields, and the 
greater the mass of the field the slower it decays 0, Ql- At a linear level, purely real modes which 
corresponds to non-damping oscillations can appear, leading to the appearance of infinitely long living 



modes, the so-called quasi-resonance modes 



53|-|55|. 



Besides that, massive scalar fields seem to behave very differently at very late times. At late times 
massless scalar and gravitational fields perturbations of Schwarzschild black holes were shown by Price 
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571 ] . to decay according to a power law of the form 

|*| ~ t- (2/+3) , (3.1) 



where I is the multipole number. Instead, massive fields have oscillatory late-time tails. One of the 
reasons, is that massive field tails appear already in Minkowsky spacetime. This is related to the fact 
that different frequencies forming a massive wave pac ket have different phase velocities. In flat spacetime 



the late-time tails of the scalar field are given by 
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(3.2) 



where [i is the mass of the field. Normally, the Minkowsky spacetime tail shows itself in the black hole 
tails at intermediate late times. In the presence of a black hole, these intermediate tails are not the final 
asymptotic behavior. In fact, at very late times, the massive scalar field decays as 
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|*| - t- 5/6 sm((it), 



(3.3) 



independently of the number I. These asymptotic tails are believed to be a resonance backscattering due 
to the curvature-induced potential. 

Here we show that in a collision of a point particle with a Schwarzschild black hole, the late-time 
tails due to the curved background, for the lowest radiatable multipoles of the massive scalar field, are 
dominant even at intermediate late times over all other contributions, namely, the quasinormal ringing 
and the Minkowsky spacetime tail. 
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3.1 Zero- frequency limit 

In the massless scalar radiation case some characteristic features of energetic collisions of point particles 
with black holes are expected, namely: the spectrum and the waveform largely depend on the lowest 
quasinormal frequency of the spacetime under consideration which works as a cutoff for the energy 
spectra; there is a non-vanishin g ze ro-frequency limit (ZFL) for the spectra, whereas for low-energy 



collisions the ZFL is zero 



49 



61 



,|6J. 
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62j, but considering 



To understand this, let us do the classic ZFL calculation for head-on collisions 
the emission of scalar radiation instead of gravitational radiation. 

The initial configuration consist of one point particle with scalar charge q p freely moving toward a 
chargeless particle, with four velocity and constant positive velocity v corresponding to a boost factor 
°f 7 = (i_ 1 ,2- l i/2 ■ At i = the two particles collide instantaneously, forming one chargeless particle and 
emitting scalar radiation. Note that the non conservation of the scalar charge is motivated by the no-hair 
theorem. Since the initial chargeless particle and the final particle do not contribute to the emitted scalar 
radiation, the stress-energy tensor for this system is simply 

T^(x, t) = & u"u"a 3 (x - vf)G(-f) , (3.4) 

7 

where the boldface denotes a three- vector. 

Using eq. (|2.5[) . the energy per solid angle and per unit frequency emitted in the wave direction 
k = k/<j is 

where the Fourier transform of the stress-energy tensor (|3.4[) is given by 

r^(k,„)= fr"""" (3.6) 

Ztt i 'y { to — v • k) 

For convenience the axes are oriented such that the motion of the initial particles is in the z axis. We 
then trivially find 

^ * «. (3.7) 



duodVt 167r 2 7 2 (1 — v cos ( 
where 9 is the angle between the particles motion direction and the wave direction. The total energy 
diverges unless a cutoff frequency ui c is introduced, being the total radiated energy given by AE = ^fjw c - 
Numerically we will see that this is indeed the case: the energy spectra for high velocities is approximately 
flat until a cutoff, which corresponds to the fundamental quasinormal frequency of the black hole. 

In order to compare the ZFL calculations against numerical results we perform a multipolar decom- 



position of the radiated energy using 



47] 



dujdQ \ ^ V doj 



where Yi m are spherical harmonics and \j are yet undetermined functions of cj. From the orthonor- 
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mality of the spherical harmonics it follows that 



dui 



d 2 E 

1 

dujd&l 



(3.9) 



Since we considered a collision along the z axis there is no dependence on the azimuthal angle </>, and 
thus only the m = modes contribute. The monopole (I = 0) contribution then reads 



dE 0C 
duj 



87T 2 7 2 V 



■ [arctan(w)]^ 



(3.10) 



In Table |3~T1 we compare the ZFL for the I = mode evaluated using Eq. (|3.10j) against ^fr-|u;=o 
obtained numerically using the point-particle method discussed in the next section. Our numerical 
results are in very good agreement with the ZFL predictions. 



E 


ZFL ana (x<7-^) 


ZFL nura (x<j- 2 ) 


1.5 


0.00938 


0.00939 


5 


0.00277 


0.00276 


10 


0.00115 


0.00113 


15 


0.000653 


0.000663 



Table 3.1: Comparison between the ZFL of the I = mode evaluated analytically and numerically, for 
E = 1.5, 5, 10, 15. The agreement is remarkable. 



3.2 Numerical Setup 

The formalism discussed in the last section is a flat-space approximation valid for the low-frequency part 
of the energy spectrum. In this section we consider a test particle with scalar charge q p and gravitational 
mass m p , and a massive scalar field coupled to matter, falling into a Schwarzschild black hole along a 
radial timelike geodesic. This is an accurate description at all frequencies in the limit where one of the 
binary components is much more massive than the other. 

In this particular case the background metric is the Schwarzschild metric which is given by 

(if 

ds 2 = -f(r)dt 2 + —— + r 2 (d9 2 + sin 2 9d(j> 2 ) , (3.11) 
f(r) 

where f(r) = 1 — ™. The test particle is described by the stress-energy tensor 

f dr 

T^= —^=q p S A (x-X(T))x fl x u , (3.12) 



where is the trajectory of the particle along the word-line parameterized by his proper time. 

Using the framework introduced in chapter [2l and specializing for the case of a test-particle radially 
infalling into a Schwarzschild black hole in four dimensions, we arrive at a wavefunction for the scalar 
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field whose evolution is given by the wave equation 



X(u,r) = f(r)S , 



(3.13) 



where the potential V is given by, 



V = f(r) fjf + 



a 1(1 + 1) 2M\ 



(3.14) 



2.121) . setting a = 0, n = and redefining the 
49(. The source term S depends uniquely on 



Note that this equation corresponds to the wave equation 
source term in a way that it coincides with the literature 
the stress-energy tensor and on the geodesic the particle follows. For massive particles, the radial timclike 
geodesies can be written as, 



dT 
dr 



E 



dT 



E 



f(r)^E 2 -l + 2M/r ' dr f(r) 



(3.15) 



where E is a conserved energy parameter. If we consider a particle with velocity Voa at infinity, then, 

E = (l_„2Jl/2 = 7- 



For a massive point particle the source S is given by (sec Appendix [Bj) . 



5 = -4^^(0,0)e^«(J)- 1 



2nr 



(3.16) 



Here, Yi m (6, </>) are the spherical harmonics and the particle velocity is given by j- = —\JE 2 — 1 + 2M/ 



The energy spectra is given by (see Chapter [5l Eq. (|5.25p ) 



dE 



— =ojV^^\X(u J ,r)\ 2 



(3.17) 



and to reconstruct the wavefunction as a function of the time t we use the inverse Fourier transform 



X(t,r) 



e- lu,t X(uj,r)duj . 



(3.18) 



To find X(oj,r) we use the Green's function technique described in chapter [2j Using Eq. (|2.24[) we 
get, at infinity, 

l(w,r^oo) = — / X th Sdr, (3.19) 

vr Jr H 

where X th is the solution of the homogeneous wave equation with the correct boundary condition at 
the horizon and W is the wronskian of the homogeneous solutions of Eq. (|3.13|) given by (|2.23l) . We 
find A{ n by solving Eq. Q3.13P with the right hand side set to zero, and using the boundary condition 
given in Eq. (|2.21[) . Matching the solution to the asymptotic solution at infinity we then find A m . For 
computational purposes, good accuracy is hard to achieve, so we keep higher-order terms of the condition 
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at infinity, using an expansion of the form: 



x r H „ e i^^r r fi + + ^ + 



(3.20) 



with, 



%M (2oj 2 - n 2 ) 
\Ju 2 - [i 2 

-iA + ifi 2 +f3 Um^Jw 2 - [i 2 2M -^cu 2 - - 2iM (cu 2 - fj, 2 ) 



2M 



C 2 



2(/3 - l)\/w 2 - At 2 - 2iM(2w 2 - /x 2 ) 
iA - 2i/3 2 + f3 (ACi yju 2 - /i 2 + 3i) - 5Ci ^ 



cj 2 — fi 2 — i 



- id (A - /3 2 + 3/3 - 2) 



A-/ 2 



2(/3 - 2) ^Juj 2 - \i 2 - 2iM (2w 2 - [i 2 ) 
4(1 - 2/3)^^-^ - 4iw 2 Ci + 4ifi 2 d 



2{P - 2)^Jw 2 - n 2 - 2iM (2oj 2 - /i 2 ) 



(3.21) 



where A = 1(1 + 1). 

When doing the numerical integration of (|3.19[) . convergence is hard to achieve. To assure convergence 
we integrate outward until a large value of r, typically r = 5000/cj. 



3.3 Numerical Results 



In Fig. 13.11 we show the energy spectra with an energy parameter of E = 1.5, for the three lowest values 
of I. As expected the ZFL predictions are in very good agreement with our numerical results and the 
fundamental quasinormal frequency for each mode act as a cutoff. 

The waveform as a function of the retarded time u = t — r* for the massless case is shown in Fig. 13.21 
At early times the wavefunction is not zero, reflecting the fact that the particle begins to fall with a non 
zero velocity. At late times the signal is dominated by the quas inormal ringing. As expected, our results 
are very similar to the ones obtained in the work of Ruffini [46( and Cardoso and Lemos 49( , where they 
computed the gravitational radiation emitted in the same process, but with E — > oo. 

The energy spectra for the massive scalar field {^M = 0.05) case is shown in Fig. 13.31 In this case the 
field mass acts as lower cutoff since no energy can be radiated for frequencies below the field mass. We can 
also see that the quasinormal frequency of the black hole acts as a upper cutoff, as in the massless case. 
Comparing these spectra with the ones obtained for the massless case, we can see that the real part of 
the quasinormal frequencies of a Schwarzschild black hole for massless and massive scalar perturbations, 
arc very similar, as pointed out in Ref. [54I ]. 

The waveforms for a scalar field of mass \iM = 0.05 at a fixed radius r* = 10M , are given in Fig. 13.41 
For the quadrupolar mode / = 2, at late times, the signal is clearly dominated by the quasinormal ringing 
with frequency uj ~ 0.48 /M. However in the lowest radiatablc modes, I = and 1 = 1, the quasinormal 
ringing does not seem to appear even at intermediate late times. Instead the signal is dominated by 
a tail of the form X ~ i~ 5 / 6 sin(/xt), independent of the angular number I. This is shown in Fig. 13.51 
where we can see that this tail fits very well the numerical curve at late times. The contribution from 
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fiM = 

1=0 

1=1 

1=2 











0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Figure 3.1: Energy spectra of the massless scalar radiation for the three lowest multipoles, for a massive 
particle falling from infinity into a Schwarzschild black hole with E = 1.5. The vertical lines correspond 
to the real part of the fundamental quasinormal mode for I = 0, I = 1 and I = 2 given respectively by, 
Mojr — 0.11, Mur = 0.29, and Mujr = 0.48. The horizontal lines correspond to the ZFL predictions. 
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Figure 3.2: Waveforms of the massless scalar radiation for the three lowest multipoles, for a massive 
particle falling from infinity into a Schwarzschild black hole with E = 1.5. Here, the wavefunction X is 
measured in units of q p . 
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the Z dependent Minkowsky tail X ~ t C+ 3 / 2 ) sin(/zi) is also shown, and it is clear that at intermediate 
late-times this contribution is not negligible. 

0.005 
0.004 



0.003 

H I J 

^ ^ 0.002 
0.001 



0.000 


Figure 3.3: Energy spectra of the massive scalar field of mass \xM = 0.05 for the three lowest multipoles, 
for a massive particle falling from infinity into a Schwarzschild black hole with E = 1.5. The vertical 
lines correspond to the real part of the fundamental quasinormal mode for 1 = 0,1 = 1 and / = 2 given 
respectively by, Muj r = 0.11, Mu R = 0.29, and Mw R = 0.48. 




It is important to point that the curvature dependent tail X ~ t 5 / 6 sin(/i£) is quite universal not 



only because it does not depend on Z, but also because it appears also in Kerr black holes 
other massive fields (Dirac 



63| . and for 



63l . |64[ and Proca [65J ) . It is thus expected that this behavior is universal for 
massive fields and does not depend on the details of the black hole horizon geometry [(35 1. Therefore, the 
signal emitted by the lowest multipoles of the massive scalar radiation does not give us much information 
about the black hole parameters, but can give us information about the field mass. This behavior was 
also found for other values of the mass /x, thus confirming our results. 
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Figure 3.5: Late-time tails of the massive scalar radiation (fiM = 0.05) for the two lowest multipoles. 
The blue line corresponds to the theoretical contribution of the Minkowsky tail X ~ £-( ; + 3 / 2 ) sin(/it), 
the purple line is the contribution of the curvature tail X ~ i~ 5 / 6 sin(/xi) and the red line is the sum of 
the two tails. 
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Chapter 4 



Superradiance in 4 + n dimensions 



si 



The phenomenon known by the name of superradiance was first proposed by Zel'dovich in 1971 [66 
He showed that a cylinder made of absorbing material, rotating about its longitudinal axis with frequency 
fi, can amplify modes of scalar radiation of frequency ui, when the condition us < mf2 (where m is the 
azimuthal quantum number) is satisfied. He then pointed out that in a Kerr black hole a similar situation 
should arise, noting that the gravitational capture of particles and waves at the event horizon replace 
the absorption occurring in the cylinder. Superradiance is thus a process where certain wave-modes can 
be amplified by scattering off the rotating black hole at the expense of the hole energy and angular 
momentum. Press and Teukolsky showed later that the maximum amplification factor, for a Kerr black 
hole, is finite for scalar, electromagnetic and gravitational radiation and it depends on the spin of the 
perturbing wave: 0.3 percent for scalar fields, 4.4 percent for electromagnetic waves and 138 percent for 
gravitational radiation 33j,[68j. Furthermore, they showed that there is no classical superradiance for 
fcrmions [(38 1. This is to be contrasted with the well known Penrose process 0, Ql- I n fact, to extract 



This is to be contrasted with the well known Penrose process [6£ 
rotational energy from the black hole, the Penrose process only needs some decay within the ergoregion, 
regardless of the spin of the particles, thus occurring also for fermions. 

-75^. For example, we can 
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33 



71 



Superradiance is responsible for many interesting effects (e.g. 
state the strong instability of compact rotating objects without event horizons but where an ergoregion 
is present [73} . Instabilities develop when the process of superradiance repeats itself ad infinitum. For 
example, a black hole can be made unstable by placing a "mirror" around it. The wave will bounce back 
and fourth between the mirror and the black hole, amplifying at each scattering. The energy extracted 
will grow exponentially until the radiation pressure destroys the mirror. This process was named black 
hole bomb by Press and Teukolsky 



33 



721 ] . One particularly interesting example of this process is the 



case where we consider a massive scalar field, with mass n, scattering off a Kerr black hole. In this case 



0H,y,y, zzi 



Another example where a 



for uj < fx, the mass fi will play the role of the mirror 
natural mirror is present is the case of asymptotically AdS spacetimes. In this case, the boundary of the 
AdS behaves as a mirror and one would expect that in some cases instabilities develop. In fact, it has 
been shown that small Kerr-AdS black holes are unstable jz8[, but not large Kerr-AdS black holes 

We will see in the next chapter that superradiance is also responsible for strong tidal effects around 
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higher-dimensional rotating black holes. Therefore, in the context of this thesis, studying the phenomenon 
of superradiance in higher-dimensional rotating black holes is an interesting topic. In this chapter we 
shall use the framework introduced in chapter [5] to compute the energy amplification of a massless scalar 
field scattering off a singly spinning Myers-Perry black hole due to superradiance. 

4.1 Superradiance 

As we have shown in chapter [2j the wave equation describing a massless scalar perturbation of a singly 
spinning Myers-Perry black hole can be decoupled into a radial, Eq. (|2.8p , and an angular equation, 
Eq. (1231) . 

Imagine now a massless scalar wave scattering off the black hole. To see what happens to the wave we 
must impose a suitable boundary condition close to the event horizon (see Appendix [A] for a discussion 
about the correct boundary condition at the horizon). Since classically nothing can escape from the black 
hole we consider only in-going waves at the horizon, 

X kjm ~Te- ikHr ', r^r Hl (4.1) 

where T is the transmission coefficient. 

At infinity we consider that the solution is given by an incident wave and a reflected out-going wave, 

X kjm ~ e-^ r * + 1le lk °° r * , t V oo , (4.2) 

where TZ is the reflection coefficient. 

It can be easily shown that the Wronskian of the solution X k j m , with the asymptotic behavior 
described above, and its complex conjugate, which is linearly independent of Xkjm, is constant everywhere 
by virtue of the field equations. Thus, using the boundary conditions (|4.1I) and (|4.2[) , we obtain after 
some algebra, 

W(r r H ) = W{r -> oo) => \K\ 2 = 1 - (1 - |T| 2 . (4.3) 

It is evident that superradiance (\7Z\ 2 > 1) will occur if, 

< uj < mQ, H . (4.4) 

4.2 Numerical Results 

Here we want to compute the superradiant amplification for various dimensions. To do so, we numerically 
integrate the homogeneous radial equation (|2.8[) , with the appropriate boundary conditions at the horizon, 
Eq. (|4.1[) . We adopt r as being the independent variable, avoiding the numerical inversion of r*(r). We 
start the integration near the horizon, = ru + er#, where e is typically 10 . Then we integrate 
outward until a large value of r = r/, typically 77 = lOOO/w. We finally find the reflection coefficient, 
Eq. (|4.3p . matching the numerical solution at r/, to the boundary condition at infinity, Eq. (|4.2j) . 
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Fig. 14.11 shows the amplification factor for several values of the modes I = m and j = 0, which are 
the modes where the superradiant amplification is most significant, in different dimensions. Unlike the 
D = 4 case, in higher dimensions, the most amplified mode is not necessarily the I = m = 1 mode. For 
example, for a = M 1 / 1+ ™, in D = 6, it is the I = m = 2 mode, and in D = 7, the I = m = 3 mode. 
This behavior is also seen in dimensions greater than seven and is confirmed by the results obtained in 
Ref. 




Figure 4.1: Amplification factor A = (1 — |7?.| 2 ) x 100% as a function of the frequency for several values 
of the modes I = to, and j = 0. 



In Table POl we show the maximum amplification factor for the mode I = m = 1, as well as the 
frequency and spin at which it occurs for D = 4, 5, 6, 7. This table shows that the peak of the maximum 
amplification factor occurs at a finite spin parameter in all dimensions. This can be seen in Fig. 14.21 
where we show the maximum amplification as a function of the spin parameter in different dimensions. To 
evaluate A max for a given a, we varied the frequency in order to find a maximum for the amplification, and 
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n 


A 









0.366 


0.380 


0.99 


1 


0.0995 


0.505 


1.21 


2 


0.00907 


0.455 


1.37 


3 


0.000596 


0.445 


1.37 



Tabic 4.1: Maximum amplification factor of the I = m = 1 mode, for D = 4,5,6,7, as well as the 
correspondent frequency and spin parameter. 



then repeated this for each value of a. The superradiant amplification depends strongly on the dimension. 
The higher the dimension, the less the wave is amplified. In D = 4 (upper- left panel), the amplification 
increases with the rotation of the black hole and approaches the maximum value, ^4 max ~ 0.366%, at 
a = 0.99A7. In D = 5 (upper-right panel) the amplification doesn't always grow as it would be naively 
expected. In fact, it grows until a ~ 1.2M 1 / 2 , and then decreases until the extremal limit a = (2M) 1 / 2 . 
This behavior is also seen in higher dimensions. In D = 6 and D = 7 (lower panels), where there is 
no upper bound on the black hole spin, the amplification factor doesn't increase without limit as we go 
to large spins. Instead, for large spins the superradiant amplification decreases and eventually becomes 
negligible. 

The behavior of the maximum amplification factor with the black hole spin, can be partially under- 
stood computing the proper volume of the ergoregion as a function of the black hole spin. The proper 
volume can be computed using 
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V = 47r / dd n / T dOi / d-d / dr g rr 9m9<i><i>T 

JO JO i=1 Jo Jr t \ i=1 



(4.5) 



where we have considered a constant time slice, we have used the reflection symmetry of the singly 
spinning Myers-Perry black hole, and we have already integrated out the <f> dependence. The integration 
limits, Ti and r/, are set noting that the ergoregion extends between the event horizon, fj = rn , and the 
ergosphere radius, r/ = r e ; i.e. , the stationary limit. This occurs where the temporal component of the 
metric (|2.2[) vanishes. Thus, we can compute r e solving the equation g t t(r e ) = 0. In Fig. 14.31 we show the 
proper volume of the ergoregion as a function of the spin parameter in different dimensions. In D = 4 
(upper-left panel) the proper volume grows monotonically with a/M , diverging for a = M. On the other 
hand, in higher dimensions, the proper volume increases with the black hole spin, for small spins, but 
eventually reaches a maximum value and then decreases monotonically with a/M x ^ l+n \ 

Comparing Fig. l4.2l and Fig. 14.31 it is evident that there exist a correlation between the proper volume 
of the ergoregion and the superradiant amplification. This explains why the maximum amplification factor 
does not grow without limit as we increase the spin. At high spins the ergoregion proper volume goes 
to zero, constraining the energy extracted from the black hole. In fact, energy extraction from the black 
hole via supcrradiance is related to the existence of an ergoregion. Inside the ergoregion negative energy 
states are possible. If one scatter a wave off the black hole, the wave can excite negative energy modes 
which will fall into the black hole and extract energy from it. If the proper volume of the ergoregion goes 
to zero, then the wave will spend less time inside the ergoregion, extracting less energy from the black 
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Figure 4.2: Maximum amplification factor A 

max & 

function of the spin parameter a for the I = m = 
l,j = mode. Top: when D=4 (left), the maximum amplification grows with the spin parameter and 
approaches A max ~ 0.366% at a = 0.99M. When D=5 (right) the maximum amplification grows with 
the spin parameter until a ~ 1.2M 1 / 2 and then decreases monotonically. Bottom: in D = 6 (left), and 
D = 7 (right) the maximum amplification factor increases monotonically until a ~ 1AM and then 

decreases for large spins. 




Figure 4.3: Proper volume of the ergorcgion as a function of the spin parameter a/M 1 ^ 1+n K 
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Figure 4.4: Angular velocity at the horizon of locally nonrotating observers £Ih = a /( r H + a2 ) as a 
function of the spin parameter a/M 1 ^ 1 " 1 "™^. 



hole and consequently, the maximum amplification factor will also asymptotically vanish. 

For a matter of completeness, we can also compare the maximum amplification with the angular 
velocity at the horizon of locally nonrotating observers £Ih = a /( r H + a 2 )- Comparing Fig. 14.21 and 
Fig. 14.41 there seems to be a better correlation between the superradiant amplification with the proper 
volume of the ergoregion than with . That is evident from the behavior of with the spin parameter 
in D = 5. Unlike the maximum amplification it grows linearly with the spin. Moreover, in D = 6 and 
D = 7 we can see that £Ih decays slower at large spins than the maximum amplification. 

An evident functional correlation between the maximum amplification and the ergoregion proper 
volume is difficult to find, and as we can see comparing Fig. 14.21 and Fig. 14.31 the maximum peak of the 
amplification factor and of the proper volume does not occur at the same value of the spin a. However, it 
is interesting to note that, for a > 1.5M 1 ^ 1+n \ ^4 max grows linearly with the ergoregion proper volume. 
In Fig. 14.51 we show the maximum amplification factor as a function of the ergoregion proper volume in 
D = 6 (left panel) and D = 7 (right panel), for large spins. Fitting the data to 

A max = ai +hV, (4.6) 

we find, for D = 6, 

oi = 0.000500012 , bi = 0.000345673 , (4.7) 
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and, for D = 7, 



a x = -0.000074633 , b x = 0.000034944 . 



(4.8) 



We checked this linear relation between the maximum amplification and the ergoregion proper for large 
spins up to D = 9, which led us to say that it should be valid for any dimension. 




10 15 
Volume/M 5 




10 12 14 
Volume/M 6 



Figure 4.5: Maximum amplification factor A max as a function of the ergoregion proper volume for a > 
1.5M 1 /( 1+ ™'. The dots are the numerical results and the straight lines correspond to the fitted curves. 
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Chapter 5 

Tides for charged interactions in 
4 + n dimensions 



In the second scenario discussed in this thesis, we consider a particle orbiting a higher-dimensional 
rotating black hole. We will see that in higher dimensions, strong tidal effects arise and should be crucial 
to determine the binary evolution. 

In fact, gravitational binaries are intrinsically complicated systems that display a wealth of interesting 
effects. One important effect, which occurs, for example, in several planet-moon systems, are tides 
generated by differential gravitational forces. Tidal forces in the Earth-Moon system have long ago locked 
the Moon in a synchronous rotation with the Earth and have increased the Earth-Moon distance 



82 



S3 



These processes, so-called tidal locking and tidal acceleration, respectively, are possible only if there is 
some dissipation mechanism in the system. Because of friction, tides extract energy from the binary, and, 
since angular momentum must be conserved, this provides a mechanism to exchange angular momentum 
between the Earth and the Moon. In the Earth-Moon case the dissipation is caused by the friction 
between the oceans and the Earth surface. In binaries containing rotating black holes the event horizon 
can play the role of a dissipatiye membrane, and tidal acceleration is known for many years under a 
different name: superradiancc 
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As stated in the last chapter, superradiance is responsible for many interesting effects, and one of them 
is the possible existence of "floating orbits" around black holes. Generically, orbiting bodies around black 
holes spiral inward as a consequence of gravitational- wave emission. When the condition for superradiance 
is met, it is possible to imagine the existence of floating orbits, i.e. , orbits in which the energy radiated to 
infinity by the body is entirely compensated by the energy extracted from the black hole 



Within 



general relativity in four dimensions, tidal effects are in gen eral completely washed out by gravitational 



wave emission and orbiting bodies always spiral inward j25| . However, when coupling to scalar fields is 
allowed, an induced dipolc moment produces a tidal acceleration (or polarization acceleration 25[), which 
might be orders of magnitude stronger than tidal quadrupolar effects. Furthermore, in theories where 
massive scalar fields are present, the coupling of the scalar field to matter can produce resonances in the 



scalar energy flux, which can lead to floating orbits outside the innermost stable circular orbit 



HQ. 
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It was recently argued via a tidal analysis framework that higher-dimensional black holes in general 
relativity should be prone to strong tidal effects [25[. One of the consequences of those studies was 



that orbiting bodies around higher-dimensional rotating black holes always spiral outward, if the tidal 
acceleration (or, equivalently, the superradiance) condition is met. 

In this chapter, we use a fully rclativistic analysis, albeit in the test-particle limit, to prove this 
behavior. For simplicity, we consider the coupling of massless scalar fields to matter around a rotating 
black hole in higher-dimensional spacetimes. We show that, for spacetime dimensions D > 5, tidal effects 



gy radiated to infinity. 
; in our case we view 



are so strong, that the energy extracted from the black hole is greater than the ener 
Higher-dimensional spacetimes are of interest in a number of theories and scenarios 
them as a proof of principle for strong tidal effects in black hole physics, without the need for resonances. 
We do not consider gravitational perturbations; gravitational effects should be subdominant with respect 
to the dipolar effects discussed here [25J]. Nevertheless, the arguments presented in Ref. 2fj] together 
with the present results show that a purely gravitational interaction also displays this phenomenon, 
which likely leads to new interesting effects in higher-dimensional black hole physics. 



5.1 "Polarization" acceleration 

The flux emitted by a particle orbiting a spinning black hole can be estimated at Newtonian level in 



terms of black hole tidal acceleration and by applying the membrane paradigm 



251. In this section, we 



generalize the computation sketched in Ref. [25| to higher dimensions and to massless scalar fields. 

Let us consider the interaction of a particle with scalar charge q p and gravitational mass m p orbiting 
a neutral central object of mass M and radius R. If the object has a dielectric constant e = e r eo, the 
particle external field induces a polarization surface charge density on the central object and a dipolc 
moment, which are given, respectively, by 



where 



(Tpoi = (3 + ri)e j3E cos>&, (5.1) 
p = n (n+3) e pR 3+n E , (5.2) 



Eo = qp 2+n , (5.3) 

"(n+3) e 0^0 



and ro is the orbital distance, Q( n +3) is the solid angle of the (n + 3)-sphere, (3 is some constant that 
depends on the relative dielectric constant of the object, and # is the polar angle with respect to the 
single axis of rotation of the central object. 

Assuming circular orbits, the tangential force on the charge q p due to the induced electric field is 
given by 

"(n+3) £ r 

Without dissipation, the dipole moment would be aligned with the particle's position vector. Here, we 
consider that dissipation introduces a small time lag r, such that the dipole moment leads the particle's 
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position vector by a constant angle <j> given by (see 25J, [83| for details) 



(5.5) 



where fin and Q are the rotational angular velocity and the orbital angular velocity, respectively. At 
first order in </>, the tangential component of the force reads 



q P p 



^(n+3) e r + "' 



(5.6) 



This exerts a torque roF$ and the change in orbital energy over one orbit reads 



1 



2 71 



2jv 



^orbital = 7— / r F$ri d-d = flr F$ 



O 4 

"(n+3) e r 



4+2n 



Q(ft-fiff)T, 



(5.7) 



where, in the last step, we used Eq. (|5.2 



5.1.1 "Polarization" acceleration of black holes 



Remarkably, the equation above qualitatively describes the energy flux across the horizon of a rotating 
black hole if one identifies with the angular velocity of the black hole and the lag r with the light- 
crossing time, t ~ R/c, where R 1+n = GdM/[(ti + 1)c 2 ], Gd is the D-dimcnsional gravitational constant, 
and M is the black hole mass j^ . Accordingly, a particle orbiting a rotating black hole in 4+n dimensions 
dissipates energy at the event horizon at a rate of roughly 



^ iTJTW "i+T 4+2» »(" ~ "gJ ■ I 5 ' 8 ) 

^(n+3)£0C *+» (n+l)*F= r o 

On the other hand, charged accelerating particles radiate to infinity according to Larmor's formula, which 



in 4 + n dimensions reads 



[this will also be derived in the next section, see Eq. (|5.48[) ] 



2„,2 



%1 

r 3-\-n 



n 4+n r% , (5.9) 



where 7 is some coupling constant. 



Tidal acceleration [25f occurs when the orbit of the particle is pushed outward due to energy dissipation 
in the central object. This is possible only if two conditions are satisfied: (i) f2 < so that Eh < and 
the energy flows out of the black hole; and (ii) \Eh\ > E^, i.e. the rate at which energy is dissipated to 
infinity must be smaller than the rate at which energy is extracted from the black hole. From Eqs. 
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and (|5.9j) . and using Q ~ r ( 3+n )/ 2 



861, we find 



1^1 /?G 



1+33 
1+ti 



7? ^ O (r,+3)(2-r») 

M T+^ (-1U-+3 



" +1 . (5.10) 



where we have assumed £7// 3> and we have defined the orbital velocity 

v = [M(n + 1)}^ fl^ . (5.11) 

At large distance, v ~ rof2, and when n = 0, we recover the standard definition, t> = (Mf2) 1//3 . 
Surprisingly, for n > 1 (D > 5) tidal acceleration dominates at large distances. This simple argument 
suggests that test particles orbiting rotating black holes in dimensions greater than five would gcncrically 
extract energy from the black hole horizon at a larger rate than the energy emitted in gravitational waves 
to infinity. As a consequence, the orbital separation will increase in time; i.e. the system will "outspiral" . 
Using the framework presented in chapter [5J in the next sections we shall prove this is indeed the case, 
by computing the linear response of a higher-dimensional spinning black hole to a test particle in circular 
orbit. 



5.2 Tidal acceleration: Analytical solution at low frequencies 

This section contains one of the main results of this thesis. We solve the wave equation, Eq. (|2.12|) . 



analytically in the low-frequency regime (see e.g. 87H9lj). obtaining a formula for the scalar flux which 



will be compared with the results obtained by a direct numerical integration of the wave equation. 

To solve the radial equation in this regime we use a well-known matching technique. We first solve 
the equation in the near-horizon regime, then we find the equivalent in the far field limit. Finally, we 
stretch and match the two solutions in an intermediate zone. This way we can construct an analytic 
expression for the radial part of the field valid in the entire spacetime. 



5.2.1 Setup 

Extending the results of the last section, we will restrict to equatorial circular orbits (i? = 0, $ = 7r/2), 
which is an unrealistic approximation in higher dimensions: generic circular orbits are unstable, with an 
instability time scale of order of the orbital period 86|. Nevertheless, our purpose here is to show that 



tidal effects can dominate, and it is not clear what the overall combined effect of tidal acceleration and 
circular geodesic motion instability is. Extending the present analysis to generic orbits and relaxing the 
test-particle approximation are interesting future developments. 

For prograde orbits around a singly spinning Myers-Perry black hole (|2.2[) the energy, angular mo- 
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mentum, and frequency of the point particle with mass m p orbiting at r = rp read [8 



3+, 



£p _ ay / (n + 1)M + r 2 - 2Mr 



7 3 + " / ^ 3 +" 1 

p r 4 y/2ay/(n + l)M + r 2 -(n + 3)Mr 



777 3(n + l) / 3 + , t l] 

p r 4 y2 av /( ? 7 + l)M + r 2 - (n + 3)Mr : 



V(« + W (r 2 - 2a^r ^ + a 2 ) 



(5.12) 



(5.13) 



J(n+1)M 

°p = — • ( 5 - 14 ) 

av /(n + 1)M + r 2 

The only nonvanishing components of the (4 + 7i)-vclocity U v of the particle on a timelike geodesic are 
given by 

, (n 2Ma 2 \ 2MaL v 
rripAr^U' = [ r 2 +a 2 + ^ fr ) E p ^ , (5.15) 

m p A r = ro U* = ^g^ + (l-4^) ip (5.16) 
As seen in chapter [5] the nonhomogeneous radial equation for the scalar field reads 



dr 2 



A 



WO = (r 2 + a2)3/2 ^ /2r ^- - ( 5 - 17 ) 



where the source term Tfe m j is given, in this case, by 

Tkmj — ~jj^Simj( n /2)Y*(ir/2, 7r/2, . . .)5(r — ro)S(mCl p — u>) , (5.18) 

which has been derived from the stress-energy tensor of the point particle in an equatorial circular geodesic 
around the black hole (see Appendix |B|) . We note that at i9 = tt/2 only higher dimensional spheroidal 
harmonics with j = are nonvanishing. This can be seen from Eq. (|2.19p . Thus, in order to calculate the 
fluxes on circular orbits, one only needs to consider terms with j = 0. In this case, the hyperspherical 
harmonics Yq are constant. 

5.2.2 Energy Fluxes 

Using the Green's function approach we can derive formulae for the energy fluxes at infinity and at the 
horizon. 

For very large values of r, Eq. (|2.24[) has the following asymptotic form: 

pik^r, poo \ In/2 



Xkjm (r -+ oo) = ^j-j- j ^ T kjm {r>)Xlf m {ri2 + a2)m dr>* = Z% m 6(w - mQ p )e^ , (5.19) 
where, using Eq. f)5 . 18|) . 



KfM S* kjm ( n/2)Y?(n/2, n/2, . . .) 

V r Q + a2r o' 



Z ?o = -a " mv ;' "' mv 1 L± ' \ ' ' q . (5.20) 
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Likewise, at the horizon we get 



X kjm (n -> -oo) = Z r «5(uj - mQ p )e- tkHr ' 



(5.21) 



where, 



z r H _ X% m (r )S* kjm (7r/2)Y?(n/2, Tr/2,...) 



/2 



The scalar energy flux at the horizon and at infinity are defined as 

« n 

E H , oo = hm / ddd<j> TT d^V^T t r , 

J i=l 

where the stress tensor reads 



(5.22) 



(5.23) 



Finally, using Eqs. (|5.19j) and (|5.21l) . we get 



(5.24) 



E 



H,oa 



kjm 



iij£l p k H ,oo\Z l 



th ,oo i2 
kjm I 



(5.25) 



The equation above shows that, il the superradiant condition kn < (w < mSlu) is met, the energy flux at 
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681. Note 



the horizon can be negative; Eh < 0, i.e. energy can be extracted from a spinning black hole 
that this condition is exactly the same as the one obtained from the Wronskian analysis, Eq. (|4.4p . In four 
dimensions, \Ejj\ "C -Eoo and the superradiant extraction is generically negligible. As we show below, in 
higher dimensions the opposite is true, \Eh\ 2> -Eoo and superradiance dominates over gravitational- wave 
emission. 



5.2.3 Solution X r k f m 
Near-horizon regime 

Let us first focus on the solution X^J m , which is regular at the horizon. 
We first make the following change of variable: 

A dh A(r) 

h = -= ? =► — = (1 - h)r = W 9 , 5.26 

r 2 + a 2 dr v ' r 2 + a 2 K ' 

where A(r) = (n + 1) + (n — l)a 2 /r 2 . Then, near the horizon r ~ rn, the radial Eq. (|2.8[) can be written 

as 

, /-, , n d 2 R . dR 

h{1 - h) ^ + il - DJl) lh + 
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P- 



A 



A(r H ) 2 h{\ - h) r 2 A(r H ) 2 (l-h) 



R = 0. 



(5.27) 



where 



P = Ld(r H + a 2 /r H ) - ma/r H , 

A = [1(1 + n + 1) + j(j + n - l)a?/r 2 H ](r% + a 2 ) 



P>* = 1 



4aV 



H 



[(n + l)r| + (n - l)a 2 ]" 

Using the redchnition = /i Q (l — h)^F(K) the above equation takes the form 

d 2 F dF 
Ml - h)—r + [c-(a + b+ l)h]— - (ab)F = , 
ah 1 an 

with 

a = a + /3 + D*-l, 6 = a + ,5 , c = 1 + 2a , 
where a and /3 must satisfy the following algebraic equations: 

-a ' ^ -0, 



£ 2 + /3(F>» - 2) + 



P 2 



A 



^(r^) 2 r 2 H A{r H y 



0, 



whose solutions read 



P 



0± 



A(r H ) ' 



±J(£», -2) 2 -4 



p2 



A 



-4(^) 2 r 2 H A{r H f 



(5.28) 
(5.29) 

(5.30) 



(5.31) 
(5.32) 

(5.33) 
(5.34) 



(5.35) 



(5.36) 



The two linearly independent solutions of Eq. (|5.31[) are F(a, 6; c; h) and h F(a + 1 — c, b+l — c; 2 — c; /i), 
where F is the hypcrgcomctric function. Convergence requires Re [c — a — b] > 0, which can be obtained 
only if the minus sign is chosen in the solutions above. In the following, we shall identify /? = /?_ and 
a = a-. The general solution of Eq. (|5.27p is then 



R(h) = Ai/i a (l - hfF(a, b, c; h) + B 1 h- a {\ - h) f) F{a + 1 + c, b + 1 - c, 2 - c; h) . 



Expanding the above result near the horizon, we get 



(5.37) 



R(h) = Aih + B^ti^^ = A ie - lkHr * + B 1 e lkHr * 



(5.38) 
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Regularity at the horizon requires B\ = 0. The near- horizon solution can be written as (see [92j) 



We can now expand this result in the low-frequency regime and for small values of a/ru, in the region 
where l-/iCl and r ^> rn, 



R 



^kjm 
r l+n/2 



T[l + 2a]r[2 - - 2(3} 



(2M)^4 T[2-D.+a- /?]r[l +«-/?] 



(5.40) 



Far-field limit 

We will follow Poisson to solve the wave equation at large distances [87[ . It is useful to rewrite the radial 
Eq. (|5.17[) in terms of the dimensionless variable z = ujt. At large distances Eq. (|5.17|) reads 



d 2 (n + l)e d Z(/ + n + l) + §(l + f) £ [i + n (^ + i 

*d^ + ~ 



z 2+n dz 



3+n 



X kjm (z) = 0, (5.41) 



where / = 1 — e/z n+1 , e = 2Mlo u+1 is a dimensionless parameter, and we used the fact that at large 
distances the eigenvalues take the form Aij m = 1(1 + n + 1). We can rewrite it in a simpler form if we 
define the quantum number J(J + 1) = /(/ + n + 1) + § (1 + § )• Solving for J, and assuming J is a 
non-negative number, we get 



In the limit e < 1, Eq. (|5.41[) reads 



J = l 



J(J + 1) 



dz 2 



Xkjm(z) = . 



(5.42) 



(5.43) 



The solution can be written in terms of a linear combination of Riccati-Bessel functions, \ f zJj+i/2( z ) 
and \fzNj + i/2{z). The requirement that X r h " m (z) be regular at the horizon demands 



X r k f m (z) = B^J J+1/2 (z), 



(5.44) 



where B is a constant. The asymptotic expansions for the Bcsscl functions are well known and read 

Bz J+1 



2 J + l / 2 T[J + 3/2] 



[l + 0(z 2 )] 



(5.45) 



At large distance, the second term within the square brackets is subdominant and we shall ignore it. 
Matching (jS^tipj) to (pT4T))) we get 



2 .7+l/2 r[J + 3/2]r[1 + 2Q]r[2 _ D# _ 2/ g] (2M) l/(2n+2) 

T[2-D. + a-0\T[l + a-0\ e (J+D/(«+V/ 2 [ + (£jJ 



(5.46) 
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The parameter A m can be extracted from the behavior of the function near z = oo. Recalling the large 
argument of the Bessel functions, X^ m (z — > oo) ~ Byj2/it sin(z — Jtt/2) , and, using Eq. (|2.21j) . it follows 
that 



^4in — 



2 J r[J + 3/2]r[l + 2a]r[2 -D t - 2/3] 
0Fr[2 - D* + a - /3]r[l + a - J] V £l/(n+1) 



J+i 



(2M) 



l/(2rc+2) 



1/2 



[1 + 0(e)] . (5.47) 



Flux at infinity 

With all of this at hand, we can now compute the flux at infinity in the low-frequency regime. From 
Eqs. (|5T2"0"j) and (|5T2"5l we get 



_ 2 i+«/2+ir[Z + n/2 + 3/2] 

2i(, 1 + l) + (n + 2)(n + 3) 



[(n + l)M] l+n/2+1 \S k3m (n/2)\ 2 \Y 3 (n/2, tt/2, . . .)p 



x r n 



(5.48) 



where we used the fact that for small frequencies (large distances) r 2 + a 2 ~ r 2 , U 1 ~ 1, and w = mf2 p 
mN /(n + l)Mro (3+n)/2 . 



5.2.4 Solution X~ m 



Let us now perform the same calculation for the solution X^ m , which satisfies outgoing- wave boundary 
conditions at infinity. The method is analogous to that already described above. However, since for this 
case the boundary condition is imposed at infinity, we do not require regularity at the horizon. In the 
limit 6<C 1, X^ m can be identified, up to a normalization constant, with 

X^ m {z) = C^H% /2 {z), (5.49) 

where Hy^^^iz) = Jj+i/2( z ) + i^j+i/2( z ) is the Hankcl function. To determine the constant C we 
match this solution in the limit z — > oo to the required boundary condition 

X^ m (z^oo)^e iz . (5.50) 

Recalling the asymptotic behavior of the Hankcl functions, we have 

C^Hfl y2 {z -> oo) ~ C S j^e i %-i) J+1 + 0(1/ z)] , (5.51) 

from which we get C = i J+1 \Jn/2. 
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The small-argument behavior of the Bessel functions reads (when J > 0) 

Jj + irf* « i) - 2 , + i/ 2 r[j + 3/2] [i + °(- 2 )] - 

AW* « 1) - - — [ zJ+1/ i l [1 + 0(z 2 )] ■ (5.52) 
At leading order and near z = 0, the function N J+1 / 2 dominates over Jj+i/2- Hence, we get 

X% m (z « 1) ~ i J ^=T[J + l/2]z- J . (5.53) 

Flux across the horizon 

We can now compute the flux across the horizon. Using Eqs. (|5.22p and (|5.25[) we get 



Z™ 

kjm 



Eh = mVlpkn 

= mk H (aq p ) 2 rl ( ) rif (2M)^ \S kjm (ir /2)\ 2 \Y 3 (tt/2, tt/2, . . .)| 2 r~^~ , (5.54) 



n l/2 

n + 1 N 



whprP T - r[i+»t/2+l/2]r[2-n»+a- f 3]r[l+ a 
wiieie i i 2r[Z+?i/2+3/2irfl+2air[2-r>.-2 



[Z+n/2+3/2]r[l+2a]r[2-Z3,-2,3] 1 

5.2.5 Ratio of the fluxes 

We can now obtain an expression for the ratio of the fluxes on the horizon and at infinity for general I, 
m, and n. Using the expressions for Eh and E^ calculated above, we find 

E H _k H r H [2 l+n/2+1 T[l + n/2 + 3/2}] 2 % ( 2 \ 2£± i i±i (^p^+d (~-ik^+i+") 



T -— (2M) WT^r 2 , (5.55) 

-Km zl+n+L \n + lj 

This can be written as a function of the orbital velocity. 

E H _ k H r H [2 l+n/2+1 T[l + n/2 + 3/2}] 2 ( 



T( x v , (5.56) 

irm 2l + n+1 1 \n + 1J K ' 

where we have used Eq. (|5.11[l at large distance. For sufficiently small orbital frequencies, such that 
the superradiancc condition is met, and the flux at the horizon is negative, we then find that the ratio 
between the fluxes grows in magnitude with ro and the particle is tidally accelerated outward. For the 
dipolar mode, I = 1, this expression is in complete agreement with the expected behavior derived from a 
Newtonian tidal analysis; cf. Equation (15. f Op . 

Note that these results were derived under the assumption of slow rotation, a -C Th- This approx- 
imation is particularly severe in the near-extremal, five-dimensional case, where th — > 0. Nevertheless, 
as we discuss in the next section, our method captures the correct scaling of the energy fluxes for any 
spin, and it even gives overall coefficients that are in very good agreement with the numerical ones in the 
slowly rotating case. This is shown in Fig. 15.11 where we compare the analytical results of this section 
with the numerical fluxes computed in the next section. 
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Figure 5.1: Comparison between the flux ratio p = E H ot /E^ t (in absolute value) calculated analytically 
and numerically, as a function of the particle velocity v for n — 0,1,2,3,4 (D = 4,5,6,7,8) and a — 
0.2M 1 /( 1+n K The straight curves correspond to the analytical formula with I — 1, and the dots are the 
numerical results discussed in section 15.31 In the slowly rotating regime, numerical results are in very 
good agreement with the analytical formula. 



5.3 Numerical Results 

The Green function approach described above can be implemented numerically u sing standard methods 
similar to the ones already discussed in the previous chapters and used in Rcf. [3,0, Hi- For given 
values of ro, a, and n, we can compute the fluxes by truncating the sum in Eq. (|5.25[) to some fc max , m ma x, 
and jmax- As discussed before, for circular orbits only j = terms give a nonvanishing contribution. 

For small and moderately large orbital velocities, the sum converges rapidly even for small truncation 
orders, and we typically set fc max = 3 and m max = 6. However, the convergence is very poor when 
the orbital velocity approaches the speed of light, i.e. when the orbit is close to the prograde null 
circular geodesic. Recall that circular orbits around Myers-Perry black holes in higher dimensions are 
unstable 



and, in particular, there is no innermost stable circular orbit for D > 4. Thus, for our 
purposes we could in principle consider particles in circular orbit up to the light ring, which exists for 
any dimension 



86j. As the particle approaches the light ring, the flux is dominated by increasingly 
higher multipoles, thus affecting the convergence properties of the series (I5.25[) . For this reason, the plots 
presented below are extended up to some value of the velocity that guarantees good convergence. 

Furthermore, for large orbital velocity and highly spinning black holes, the zeroth order angular 
cigenfunctions and the corresponding eigenfrequencics (|2.14[) might not be accurate. Therefore, when 
au > 1, we have used exact numerical values of Akj m obtained by solving Eq. (|2.9|) with the continued 



fraction method described in chapter [5] and in Ref. 4l|. We note, however, that Eq. (|2.14p reproduces 
the exact results surprisingly well, even when au> ~ 1. 

We checked our method by reproducing the results of Ref. [3| for the massless case in four dimensions. 
In addition, we can compute the energy flux in any number of dimensions. The fluxes Eh and Eoq for 
D = 5 and D = 6 are shown in Tables 1531 and I5~2l for ?'o = lOrjj. We show the total flux as well as the 
first multipolar contributions. 
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k 


m 


3 


ro/r H 


E H {aq p ) 2 


Eoo{aq p ) 


\Eh\/ Eoo 







1 





10 


-1.4759 x 10" 


-y 


1.5790 x 10" a 


0.9347 







2 





10 


-8.3175 x 10" 


n 


1.6288 x 10~ 10 


0.5107 







3 





10 


-3.1691 x 10" 


12 


1.0080 x 10- 11 


0.3144 


1 




1 





10 


-1.2718 x 10" 


14 


5.1192 x 10" 16 


24.844 




10 


-3.1248 x 10" 


-9 


3.5052 x 10~ 9 


0.8915 



Tabic 5.1: Fluxes across the horizon and to infinity for n = 1 [D = 5), a — M 1 ^ 1+n \ and tq = 10r#. In 
the last row we show the total flux obtained summing up to fc max = 3 and m max = 6. 
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ro/r H 


E H (aq p ) 2 


Eoo(aq p )~ 


2 


\Eh\/ Eoo 







1 





10 


-4.1768 x 10- 


12 


1.3915 x 10 


-14 


300.155 







2 





10 


-3.0159 x 10- 


13 


3.8070 x 10 


-16 


792.180 







3 





10 


-1.3677 x 10- 


14 


4.8651 x 10 


-18 


2811.35 


1 




1 





10 


-2.8126 x 10- 


16 


1.0497 x 10 


-22 


2.6796 x 10 6 




10 


-8.9858 x 10" 


12 


2.8601 x 10 


-14 


314.176 



Table 5.2: Same as in Table Q but for n = 2 (D = 6). 



Tables 15.11 and 15.21 confirm our analytical expectations that the behavior for n > 1 (D > 5) is 
qualitatively different: the energy flux across the horizon is larger (in modulus) than the flux at infinity. 
This is shown in Fig. 15.21 where we compare the flux ratio p = Ejf 1 /E™ as a function of the orbital 
velocity v for a = 0.99M 1 ^ 1+n ^ ) in various dimensions. This figure is analogous to Fig. 15.11 but for 
a = 0.99A/ 1 /( 1+ ™', i.e. a regime that is not well described by the analytical formula (|5.56[) . For D = 4, we 
find the usual behavior; i.e. the flux at the horizon is usually negligible with respect to that at infinity, and 
the ratio decreases rapidly at large distance. The case D = 5 marks a transition, because p is constant 
at large distance. This is better shown in the left panel of Fig. 15.31 On the other hand, for any D > 5 
the flux across the horizon generically dominates over the flux at infinity. 
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Figure 5.2: The flux ratio p = E^/E^ (in absolute value) as a function of the particle velocity v 
defined in Eq. ([STTT1) for n = 0, 1, 2, 3, 4 (D = 4, 5, 6, 7, 8) and a = 0.99M 1 /( 1 +«). 



In Fig. 15.31 we show the flux ratio p for some selected value of the spin parameter a in five dimensions 
(left panel) and in six dimensions (right panel). When D = 5, the ratio is constant in the small v region, 
and it approaches unity in the extremal limit, a — > \JlM. As shown in the right panel of Fig. 15.31 
when D = 6 there exist some orbital velocity for which — p = 1, corresponding to a total vanishing flux, 
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1.0 -j 




Figure 5.3: The ratio p = Ejj ot / E^ as a function of the orbital velocity defined in Eq. (|5.11[) for several 
values of a. Left panel: when D = 5, the ratio is constant in the small v region, and it approaches unity 
in the extremal limit, a — > \J2M . Right panel: when D = 6, the flux at the horizon can exceed the flux at 
infinity. For each curve, the intersection with the horizontal line corresponds to a floating orbit, — p = 1. 
Note that, at large orbital velocity, the superradiant condition is not met and Eh > 0. 



Eh + Eoo = 0. These orbital frequencies correspond to "floating" orbits [33, [74[. Although in the right 
panel of Fig. 15.31 this is shown only for a/M 1 / 3 = 0.1,0.2,0.3, we expect this to be a generic feature 
also for larger values of the spin. The poor convergence properties of the series (|5.25[) prevent us from 
extending the curves to larger values of v, where floating orbits for a > 0.3M 1 / 3 are expected to occur. At 
smaller velocity, the energy flux contribution dominates, and the motion of the test particle is generically 
dominated by tidal acceleration. Similar results can be obtained for any D > 6. 
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Chapter 6 

Conclusions 



The main motivation of this thesis was to study what kind of effects and signatures can arise in black 
hole spacetimcs due to the perturbation of a test particle coupled to a scalar field. The interaction of 
black holes with fields is an interesting topic by itself and also fundamental for the understanding of 
astrophysical processes, and scenarios such as the AdS /CFT duality and TeV-scale gravity. 

Through the use of black hole perturbation theory, we computed the massive scalar radiation emitted 
by a test-particle falling radially into a Schwarzschild black hole. The energy spectra suffers a cutoff at the 
massive scalar held fundamental quasinormal frequency, thus showing the relevance of the quasinormal 
modes of a black hole. The signal for massless and massive scalar radiation was shown to be quite different. 
As for the massless case, at intermediate late times, the lowest multipoles of the massive radiation were 
expected to be dominated by a quasinormal ringing. Instead, the signal was shown to be dominated by 
an oscillatory tail, which decays slower than any power law. Due to the power-law decay of massless 
fields at very late-times, this behavior implies that massive perturbations decay slower than any massless 
perturbation. Furthermore, no information about the black hole parameters can be extracted from the 
lowest multipoles signal of massive scalar radiation in this kind of processes, since these oscillatory tails 
are quite universal and do not depend on the black hole geometry. 

In the second part of this thesis, we computed the amplification of massless scalar waves scattering 
off a singly spinning Myers-Perry black hole due to superradiance. This phenomenon also arise in higher 
dimensions and in dimensions greater than five, where there is no upper limit on the black hole spin, the 
amplification factor does not increase without limit when we increase the spin parameter. It was then 
conjectured that this behavior is due to the asymptotic behaviour of the ergoregion proper volume which, 
as for the amplification factor, goes to zero for very large spins. 

Finally, we concluded with the main result of this thesis, by computing the rate at which the energy 
is extracted from a singly spinning, higher-dimensional black hole when a massless scalar field is coupled 
to a test particle in circular orbit. We showed that, for dimensions greater than five and small orbital 
velocities, the energy flux radiated to infinity becomes negligible compared to the energy extracted from 
the black hole via superradiance. 

Although we considered scalar-wave emission, we expect our results to be generic in higher dimensions. 
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In particular, superradiance should be a dominant effect also for gravitational radiation. At leading order, 
the ratio \Eh\/Eoo for gravitational radiation should scale with the velocity as described by Eq. (|5.56|) . 
The dominant quadrupolc term (I = 2) reads [25 1 

\Ejj\ (r»-l)(r» + 5) 

— : r*j V " + 1 

By comparing the formula above to Eq. (|5.56|) with 1 = 1, we note that dipolar effects are dominant over 
their quadrupolar counterpart. Nevertheless, even in the purely gravitational case, tidal acceleration and 
floating orbits around spinning black holes are generic and distinctive effects of higher dimensions. 

In principle, gravitational waveforms would carry a clear signature of floating orbits {t^. [7^. Docs 
floating or these strong tidal effects have any significance in higher-dimensional black hole physics? We 
should start by stressing that circular geodesies in higher dimensions are unstable, on a time scale 



comparable to the one discussed here 



86j; however, our analysis suggests that, while more pronounced 



for circular orbits, tidal acceleration is generic and in no way dependent on the stability of the orbit under 
consideration. Wc are thus led to conjecture that tidal effects are crucial to determine binary evolution 
in higher dimensions. It is possible that tidal effects already play a role in the numerical simulations of 



the kind recently reported in Refs. [6|, [93|, [94( , but further study is necessary. One of the consequences of 



our results for those types of simulations is, for instance, that in higher-dimensional black hole collisions 
the amount of gravitational radiation accretion might play an important role. It would certainly be an 
interesting topic for further study to understand tidal effects for generic orbits, and to include hnite-size 
effects in the calculations. 

Such effects in higher dimensional spacetimes were discussed in this thesis for the first time. In the 
formation of micro black holes at the LHC, if ever recorded, strong tidal effects should certainly play 
a relevant role. This could also be the case in the AdS /CFT correspondence. Does this kind of tidal 
effects exist in AdS backgrounds in D + 1 dimensions? What other phenomenon could arise do to these 
strong tidal effects and the presence of massive fields? The understanding of the physics involved when 
dealing with black holes in four and higher dimensions is far from over and further work in this direction 
is mandatory. 
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Appendix A 

Regularity of the boundary 
condition at the horizon 



Physically motivated boundary conditions are imposed by requiring that using well-behaved coordinates 
at the horizon, i. e. , coordinates that remove the singularity at the horizon, we have neither singular nor 
identically zero fields. 

First we note that the Boyer-Lindquist coordinates are singular at the horizon. It takes an infinite 
coordinate time for any particle to fall into the black hole, t — > oo as r — > r#, and the world lines are 
infinitely twisted around the horizon, <f> — > oo as r — > m- We can remove the coordinate singularity 
transforming the Boyer-Lindquist coordinates to the so-called Kerr "ingoing" coordinates Q. In place 
of t and 4> we introduce the new variables, 

av = at H — — ar, 

d(f> = d<f> + ^dr. (A.l) 

By a direct substitution of the relations (|A.1[) into the line element (|2.2[) . the Myers-Perry metric in these 
new coordinates can be brought to the form 

, , A -a 2 sin 2 § , 2 „, , 2a(r 2 + a 2 - A) sin 2 i? , ~ (r 2 + a 2 ) 2 - Aa 2 sin 2 ?? . 2 „ , 72 
ds 2 = dv 2 + 2dvdr dvdcf) + '-— sm 2 fide/) 2 

Zj 2-i 

- 2a sin 2 Mrd^> + Y,dd 2 + r 2 cos 2 Mfl 2 n . (A.2) 

This metric is clearly non-singular at the horizon, but it is singular at S = 0, or equivalently at r = 
and <j) = 7r/2, which is the only true singularity of the Kerr-Myers-Perry space-time. 

Under these new coordinates the scalar field equation is separable using the ansatz 

tp(v, r,#,4>)=J2 [ dwe m *-^'R{r)S lmj (0)Yj . (A.3) 



l,m,j ' 
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Inserting the Kerr coordinates in the form, 



v = t+ I — — — dr 



^dr, (A.4) 



into (|A.3[) . the relation between R(r) and R(r) is easily obtained, 



l..m,j l.m,j ' 



=>fl(r) = i?(r)e l / ( " 2+ ° a""" dr . (A.5) 

At the horizon we can have either ingoing our outgoing waves: 

R~e ±ikHT % (A.6) 

where kn = w — mfifj. Requiring the regularity of R(r) around r = Th lead us to the correct boundary 
condition: 

R(r -> r H ) ~ e ±^Hr, e ij s dr _ gl j s ^ (A. 7) 

where we have used the definition of the tortoise coordinate, dr/dr* = A/(r 2 + a 2 ). Since A — > as 
r —t rn , we must have ± (r 2 + a 2 ) fc# + (r 2 + a 2 ) fc# — ^ as r — ^ rn , to assure the regularity of the solution. 
The regular boundary condition is thus, 

j?~ e ~ ikHr *. (A.8) 

Another way to see that (|4.1|) represent ingoing waves at the horizon, is checking the radial group 
velocity of the wave. For a well-behaved observer at the horizon the group velocity must be negative; i.e. 
he must see the wave progressing toward the black hole. 

The group and phase velocity of the regular solution (|A.8|) . close to the horizon, are given, respectively, 

by 

duo , . s 

-1, (A.9) 



dkn 



V phaS e=~^- = -) • (A.10) 



H 



We see that the group velocity (jA.9[) is negative, which means that for any local observer the waves are 
progressing towards the horizon. However, for an observer at infinity, Vph ase is the relevant velocity. We 
see by (|A.10[) that if m ^ H > 1 (the superradiant frequencies), the wave will appear as emerging from the 
black hole. For the observer at infinity the wave is amplified after scattering off the black hole. 
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Appendix B 



Derivation of the source terms TJ m j 
and S 



Source T) 



Imj 



The trace of the stress-energy tensor for a point particle in an equatorial circular orbit can be written in 
the form 

n 

T = -^k^S(r - r )5(tf - n/2)5(0 - Q p t) [J 5(9, - tt/2) , (B.l) 



i=i 

where, using the line element (|2.2[) , the square root of the metric determinant takes the form 

n-l 

V^9 = Ssin^r n cos"i9]^[sin6i l "- 1 . (B.2) 



To separate the wave equation, Eq. (|2.5[) . we define 

<p(t, r,0, 0) = Yl J duje im ' t> - iut R(r)S lTnj (d)Y j , (B.3) 



l,m,j 

and 



aST=^ J dwe im ^- iuJt T lmj Si mj (-d)Yj . (B.4) 

l,m,j 



Upon direct substitution of Eq. (|B.3[) . and Eq. (|B.4|) into the wave equation, we get the separate 
angular and radial equation, Eq. (|2.8[) and Eq. (|2.9D . The source term T; mj can then be obtained 
inverting Eq. (|B.4|) . 



Using (|B.1[) and (|B.2|1 , we can write the function ST in the following form 

n 

ST = - — Qp , S(r -ro)6(6- n/2)6(4> - n p t) J] 5{6 i - tt/2) . (B.5) 

£/*r»smtf cos"t9n™=i sm0" 1 fj[ 
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Then, using (|B.4|) . we get 



q p a 



-i . on _i S(r - r )S(d - tt/2)5(0 - Q p t) J] 6(9, - tt/2) 



i=l 



?7* r™ sin $ cos™ $ IX™=i sm 
= ]T J djuij m *- iut T lmj S lmj (V)Yj- 

l.m,j 

Multiplying both sides by S*, m ,j,Y*,e~ lm & and integrating on the [n + 2)-sphere, we get 
q P aS?, m , .,(tt/2)1;; (tt/2, tt/2, . . .)e"™'<V 



(B.6) 



Z7* r r ' 



-6(r - r ) = / dwe-^j 



(B.7) 



Note that, for real frequencies, the spin-weighted harmonics are orthogonal on the sphere. We chose the 
normalization, J dfle q ^ m ~ m 'Si m jS*, m ,j/YjYp = 5i>i5 m i m 5jij : which means that 



/n 
tMsintf cos"tf JJ d6 { wa-fff* Si mj S^ mj YjY* = 1 
i=l 



(B.8) 



Finally, multiplying both sides by e lcJ * and integrating in the time variable, we get the expression for 



Ti mj = --^St mj (n/2)Y*(n/2,Tr/2, . . .)5(r r )S(mQ p u>) . 



jjt r 



(B.9) 



Source S 

For the case of a radial geodesic in a Schwarzschild geometry, the trace of the stress-energy tensor for a 
point particle can be written in the form 



T 



-gu* 



5(r - R(t))6(6 - 0)6((f) - 0) 



where the square root of the metric determinant takes the form 



(B.10) 



-g = r sin I 



(B.H) 



The nonhomogeneous equation for the scalar field can written in the form, 



X lm {r*) = ^-T lm = f(r)S, 
r 



where f(r) = 1 - 2M/r and S = T hn /r. 

Once more, to separate the wave equation we define 



l T=-^=Y] ( duje- lult T lm Y lm (6, 
\2tt . J 



(B.12) 



(B.13) 



The source term TJ m can then be obtained inverting Eq. (|B. 13|) . 
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Using (|B.10|) and (|B.11[) . we can write the function r 2 T in the following form 



r 2 T 



U* sin (9 



5{r - R(t))5(9 - 0)5(<f> - 0) 



Then, using (|BTT3|) . we get 



U* sm6 



5(r - R(t))5(6 - 0)S(4> - 0) 



J duje- lult T lm Y lm (e it 

l,m 



Multiplying both sides by Y^ m , and integrating on the sphere, we get 



-5(r - R(t)) = / due^Ti 



Finally, multiplying both sides by e tLJ * and integrating in the time variable, we 



Tu 



J=^(0,0) J dt{^)-'e^S(r-R(r)) = 

J=5^(0,0) J dr{^-){f T r'e^5{r~R{r)) 

-^ = ^(0,0)e^)(^)- 1 . 
2n dr 



So, the source term S can be written, 



5 
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